home *** CD-ROM | disk | FTP | other *** search
/ Die Ultimative Software-P…i Collection 1996 & 1997 / Die Ultimative Software-Pakete CD-ROM fur Atari Collection 1996 & 1997.iso / g / gnu_c / gpplib22.zoo / libsrc / xintdbl.cc < prev    next >
Encoding:
C/C++ Source or Header  |  1993-11-30  |  3.1 KB  |  143 lines
  1. /* 
  2. Copyright (C) 1988, 1993 Free Software Foundation
  3.     written by Doug Lea (dl@rocky.oswego.edu)
  4.  
  5. This file is part of the GNU C++ Library.  This library is free
  6. software; you can redistribute it and/or modify it under the terms of
  7. the GNU Library General Public License as published by the Free
  8. Software Foundation; either version 2 of the License, or (at your
  9. option) any later version.  This library is distributed in the hope
  10. that it will be useful, but WITHOUT ANY WARRANTY; without even the
  11. implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  12. PURPOSE.  See the GNU Library General Public License for more details.
  13. You should have received a copy of the GNU Library General Public
  14. License along with this library; if not, write to the Free Software
  15. Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
  16. */
  17.  
  18. // Routines for converting between Integers and doubles.
  19. // Split up into a separate file to avoid Integer.o's need
  20. // for libm.a on some systems (including SunOS 4).
  21.  
  22. #include <xinteger.h>
  23. #include "xinteger.hp"
  24. #include <float.h>
  25. #include <math.h>
  26. #include <limits.h>
  27.  
  28. #ifndef HUGE_VAL
  29. #ifdef HUGE
  30. #define HUGE_VAL HUGE
  31. #else
  32. #define HUGE_VAL DBL_MAX
  33. #endif
  34. #endif
  35.  
  36. // convert to a double 
  37.  
  38. double Itodouble(const IntRep* rep)
  39. { 
  40.   double d = 0.0;
  41.   double bound = DBL_MAX / 2.0;
  42.   for (int i = rep->len - 1; i >= 0; --i)
  43.   {
  44.     unsigned short a = I_RADIX >> 1;
  45.     while (a != 0)
  46.     {
  47.       if (d >= bound)
  48.         return (rep->sgn == I_NEGATIVE) ? -HUGE_VAL : HUGE_VAL;
  49.       d *= 2.0;
  50.       if (rep->s[i] & a)
  51.         d += 1.0;
  52.       a >>= 1;
  53.     }
  54.   }
  55.   if (rep->sgn == I_NEGATIVE)
  56.     return -d;
  57.   else
  58.     return d;
  59. }
  60.  
  61. // see whether op double() will work-
  62. // have to actually try it in order to find out
  63. // since otherwise might trigger fp exception
  64.  
  65. int Iisdouble(const IntRep* rep)
  66. {
  67.   double d = 0.0;
  68.   double bound = DBL_MAX / 2.0;
  69.   for (int i = rep->len - 1; i >= 0; --i)
  70.   {
  71.     unsigned short a = I_RADIX >> 1;
  72.     while (a != 0)
  73.     {
  74.       if (d > bound || (d == bound && (i > 0 || (rep->s[i] & a))))
  75.         return 0;
  76.       d *= 2.0;
  77.       if (rep->s[i] & a)
  78.         d += 1.0;
  79.       a >>= 1;
  80.     }
  81.   }
  82.   return 1;
  83. }
  84.  
  85. // real division of num / den
  86.  
  87. double ratio(const Integer& num, const Integer& den)
  88. {
  89.   Integer q, r;
  90.   divide(num, den, q, r);
  91.   double d1 = q.as_double();
  92.  
  93.   if (d1 >= DBL_MAX || d1 <= -DBL_MAX || sign(r) == 0)
  94.     return d1;
  95.   else      // use as much precision as available for fractional part
  96.   {
  97.     double  d2 = 0.0;
  98.     double  d3 = 0.0; 
  99.     int cont = 1;
  100.     for (int i = den.rep->len - 1; i >= 0 && cont; --i)
  101.     {
  102.       unsigned short a = I_RADIX >> 1;
  103.       while (a != 0)
  104.       {
  105.         if (d2 + 1.0 == d2) // out of precision when we get here
  106.         {
  107.           cont = 0;
  108.           break;
  109.         }
  110.  
  111.         d2 *= 2.0;
  112.         if (den.rep->s[i] & a)
  113.           d2 += 1.0;
  114.  
  115.         if (i < r.rep->len)
  116.         {
  117.           d3 *= 2.0;
  118.           if (r.rep->s[i] & a)
  119.             d3 += 1.0;
  120.         }
  121.  
  122.         a >>= 1;
  123.       }
  124.     }
  125.  
  126.     if (sign(r) < 0)
  127.       d3 = -d3;
  128.     return d1 + d3 / d2;
  129.   }
  130. }
  131.  
  132. double
  133. Integer::as_double () const
  134. {
  135.   return Itodouble (rep);
  136. }
  137.  
  138. int
  139. Integer::fits_in_double () const
  140. {
  141.   return Iisdouble(rep);
  142. }
  143.